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We derive explicit expressions for the multi-detector ^-statistic metric applied to short-duration 
non-precessing inspiral signals. This is required for template bank production associated with co- 
herent searches for short-duration non-precessing inspiral signals in gravitational-wave data from 
a network of detectors. We compare the metric's performance with explicit overlap calculations 
for all relevant dimensions of parameter space and find the metric accurately predicts the loss of 
detection statistic above overlaps of 95%. We also show the effect that neglecting the variations of 
the detector response functions has on the metric. 



I. INTRODUCTION 

Inspiral signals are thought to be the most promising 
source of gravitational-waves (GWs) for second genera- 
tion GW detectors. Depending on the rate of merger 
events, the Advanced LIGO, Advanced Virgo GW detec- 
tor network operating at design sensitivity will be able 
to detect between 0.4 and 400 binary neutron star coa- 
lescences per year [1]. Underlying these numbers there 
is an assumed threshold on the network signal-to-noise 
ratio (SNR) at which a signal is "detectable" (i.e., has 
a false alarm probability below some established value). 
It has been shown that, among the different matched- 
filter based search strategies, coherent templated searches 
for these signals can reduce the false alarm rate for the 
same network SNR compared to coincident templated 
searches [2-5]. Thus it is attractive to prepare coherent 
searches for when the advanced detectors come online in 
order to maximize the number of detected events. 

The .F-statistic was originally derived as a single de- 
tector detection statistic associated with searching GW 
data for signals from rotating neutron-stars [6], and was 
extended to multiple-detector analysis in [ ]. However, it 
is equally applicable to coherent searches for GW signals 
from inspiralling compact objects [2-5], due to the phys- 
ical similarity of the two emitting systems. The signals 
from both types of systems can be modelled as GW emis- 
sion from a rotating quadrupole moment. Both signals 
can be characterized by four extrinsic parameters that af- 
fect the amplitude, polarization, and phase offset of the 
waveform, an extrinsic parameter that sets a reference 
time for the signal, and intrinsic parameters that affect 
the phase and amplitude evolution of the waveform. 

In performing templated matched-filter searches for 
GWs, one is always faced with the question "what tem- 
plate waveforms should the data be filtered against?" 
With regards to searches for inspiral signals in single 
detector GW data, this question has been investigated 
within a geometric formalism. Specifically, a distance 
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measure can be defined on the parameter space based on 
the "mismatch" between waveforms from different pa- 
rameter space points [7]. This was initially derived for 
the two dimensional mass space for stationary phase ap- 
proximation (SPA) inspiral waveforms expanded to New- 
tonian order in the amplitude and 1.0 post-Newtonian 
(PN) order in the phase, where the effects of the ob- 
jects' spins were neglected. This has been extended to 
3.5 PN order for the "non-spinning" contributions to the 
phase [8, 9]. In addition, a higher dimensional metric has 
been obtained that includes the "spin" contributions to 
the phase, up to 2.0 PN order, for the case where the 
objects' spins are aligned with the orbital angular mo- 
mentum [10]. 

There have been several pieces of work that have been 
closely related to deriving the multi-detector /"-statistic 
metric for short-duration non-precessing inspiral signals. 
The first was the derivation of the mismatch metric for 
coherent searches of short-duration non-precessing inspi- 
ral signals based purely on the Newtonian order inspiral 
phase model [11] and built on the formalism of [2, 7], 
which was later extended to cover the phase expanded at 
2.5 PN order [12]. Another was the derivation of the 
multi-detector /-statistic metric for rotating neutron- 
stars [ ]. In addition, there was the computation of 
the Fisher matrix for the network SNR of known and un- 
known waveforms of short- and long-duration, focusing 
on obtaining explicit expressions for the angular resolu- 
tion of a GW detector network [14]. Finally, the most 
closely related work showed parameter recovery accura- 
cies based on the Fisher matrix applied to inspiral and 
inspiral-merger-ringdown waveforms observed by detec- 
tor networks [15], although the derivation of the Fisher 
matrix was was not presented. There has been no equiva- 
lent published derivation of the multi-detector .F-statistic 
metric for short-duration non-precessing inspiral signals 
including both the amplitude model, the phase model, 
and the directional derivatives effects of detector re- 
sponses. This is what we derive here to 3.5 PN order 
in the inspiral phase. This metric is required for deter- 
mining how to arrange templates that would cover the 
four dimensional sky-location and mass space of a coher- 
ent search. 
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Previous coherent searches for short- duration non- 
precessing inspiral signals have been based on one of three 
methods. They have either relied on the sky position to 
be known precisely [16] or known to some degree and 
tiled by detector triangulation arguments [17]. These 
have both done a templated search on the mass parame- 
ter space in an ad hoc way based on mass space coverings 
associated with a single detector [ ]. A third approach 
has been hierarchical [18], relying on coincident searches 
of single detector data with their associated mass space 
coverings to decide what points in the mass space are 
followed-up coherently. The metric derived here could 
be used as the starting point for determining separately 
a template covering of the sky as well as the mass space 
covering for a template bank associated with a coherent 
search. 

Following the formalism laid out for computing the 
multi-detector F-statistic in Refs. [13, 19], this work is 
organized as follows, Sec. II identifies the form of the GW 
signal from rotating non-precessing quadrupole moments 
as seen in a GW network, Sec. Ill summarizes the formu- 
lation of the multi-detector F-statistic, Sec. IV outlines 
the approximations appropriate when applied to short- 
duration (i.e., much less than one day) non-precessing 
inspiral signals, Sec. V derives the metric for the co- 
herent multi-detector F-statistic for short-duration non- 
precessing inspiral signals, and Sec. VI shows tests of this 
metric. 



II. OBSERVED GW SIGNAL FROM ROTATING 
NON-PRECESSING QUADRUPOLE MOMENTS 

To start with, let us identify the parameters that will 
affect how a generic GW signal from a rotating, non- 
precessing quadrupole moment is observed by a GW 
detector. These parameters can be separated into two 
classes, intrinsic parameters, which affect the time evo- 
lution of the waveform and we will elaborate further on in 
Sec. IV, and extrinsic parameters, which affect the polar- 
ization, amplitude, and the phase and time offsets. The 
extrinsic parameters can be further subdivided into two 
classes, those that can be measured analytically within 
the matched-filtering process, and those that must be 
searched over by separate filters. As we will see, the ex- 
trinsic parameters that can be measured analytically are 
the extrinsic amplitude ho, the inclination angle i be- 
tween the line of sight and the total angular momentum 
of the emitting system, the reference phase 0o, and the 
polarization angle ip, which is a rotation between the ra- 
diation frame of the GW and the frame of the detector 
about the direction of propagation —h. 

With those definitions of the extrinsic parameters, we 
give a signal model that describes how a generic GW sig- 
nal from a rotating, non-precessing quadrupole moment 
will be observed in detector Y . A generic propagating 
GW signal can be described in terms of two polarizations 



in general relativity, 

h:=h + -ih x , (1) 

where the h+ and h x waveforms are out of phase by 90°. 
Thus, these waveforms can be written in terms of the 
intrinsic waveforms h c (t) and h s (t) as 

h = saf+h c (t)-is>/ x h s (t). (2) 

The intrinsic waveforms can be further decomposed into 
an amplitude piece A(i) and phase piece 0(i), 

/i c (f) := A(t)cos[0(t) + o ], 
h s (t) := A(t)sm[0(i)+0 o ], 

where A(t) and <j>(t) will depend on the details of the 
emitting system. The polarization amplitudes associated 
with the different polarization waveforms are functions 
of the extrinsic amplitude h and the inclination angle t, 

&f+ := + cos 2 t), s/ x ho cos l. (4) 

The waveform as seen by detector Vcan be obtained by 
taking the real part of this complex waveform projected 
onto the complex detector response F Y := FT + iF Y , 

s Y (t)=${hF Y e- a ^), (5) 

where we give explicit expressions for F Y and F Y later 
in this section. Expanding h, F, and the phase terms of 
cosine and sine waveforms of (5), we find 

s Y (t) = 

cos <f>o cos 2-0 — sd x sin 0o sin 2ip) F Y (t)h Y (t) 

+ cos 0o sin 2-0 + stf x sin0 o cos 2^) F Y (t)h Y {t) 

+ sin 0o cos 2-0 — srf x cos 0o sin 2^) F Y (t)h Y (t) 

+ (—^4- sin 0o sin 2^ + stf x cos0o cos 2-0) F Y (t)h Y (t). 

(6) 

This can be separated into a sum over four 
detector-independent amplitude parameters {A 11 } and 
four detector-dependent polarization-weighted waveforms 

Kit)}, 
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s Y (t) = J2^h Y (t). (7) 

It is readily apparent that the amplitude parameters are 
defined as 

A 1 := cos 0o cos 2ip — ,e/ x sin 0o sin 2-0, 

A 2 := srf-y cos 0o sin 2ip + s^ x sin 0o cos 2%p, 

, (8) 
A :=—£/ + s'm(j)oCOs2ip — £!/ x cos(t)osm2il;, 

A A := — sin 0o sin 2ip + g/ x cos O cos 2tp, 
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while the polarization-weighted waveforms are defined as 



III. THE JF-STATISTIC 



hl(t) 


:=F Y (t)h c (t- 




h Y 2 (t) 


:= F Y (t)h c (t - 


t Y ), 


h{(t) 


:=F Y (t)h s (t- 




hj(t) 


:=F Y (t)h s (t- 





(9) 



Turning our attention to the detector polarization re- 
sponses, these characterize the response of an arbitrary 
GW detector for signals that satisfy the long wavelength 
limit approximation [ ]. They can be defined as the 
double contraction of two tensors [19], 



F Y (t) :=e) 



(10) 



where dL(t) is the detector response tensor and {e+ x } 
are the polarization-independent basis tensors of the radi- 
ation frame. For an interferometric detector, the detector 
response tensor is given by 

dUt) = \ {*?(*) ® *Tw - %m ® *!(*)}.. • (ii) 

Here, V[ is the unit vector pointing along interferome- 
ter Y's first arm away from the interferometer's vertex. 
Similarly, 1% is the unit vector pointing along interfer- 
ometer y's second arm away from the interferometer's 
vertex. The polarization-independent basis tensors are 
defined as 



£ + : = 

Jj 



1 1 (8 i - fj (8 77 1 



(12) 



given in the radiation frame {£, 7}, —n}, where — h is the 
direction of propagation, and {£, 77} are basis vectors in 
the wave-plane (i.e., the plane perpendicular to direction 
of propagation) . The basis vectors £ and 77 can be defined 
with respect to h as 



; n x z f 
n x z\ 



(13) 



In a fixed reference frame centered at the geocenter, 
where 

n — (cos 5 cos a, cos (5 sin a, sin (5), (14) 

the wave-plane basis vectors are 

£ = (sin a, — cos a, 0), (15) 
77 = (— sin (5 cos a, — sin (5 sin a, cos 5). (16) 

We have now defined all of the quantities that are used 
to convert a GW signal from an arbitrary non-precessing 
rotating quadrupole moment source to the signal seen by 
a GW detector. The remaining details of the signal will 
depend on the specifics of the emitting system. 



The likelihood ratio of a signal s being in the data of 
a network of GW detectors x is given as 



A(x; s) 



P(x\s) 
P(x\0) 



exp 



(x\s) - -(s\s) 



(17) 



where (a\b) := XV( flY I" ) anc ^ the definition of the 
noise- weighted inner product (a Y \b Y ) depends on the de- 
tails of the waveform being studied. We define this for 
inspiral signals in (26) of Sec. IV. Using the signal model 
from (7), (17) can be written as 

lnAfo s) = A"x„ - ^A"M^A U , (18) 



where x L 



(x\h^) and M.^ := {hfj\h v ). In matrix 



form, M^u is block diagonal, 



May 



( A C 

C B 

A C 

\ C B 



(19) 



jl V 



due to the orthogonality of the sine and cosine intrinsic 
waveforms. Here, 



A:=Y,Fl(t)Fl(t)(h Y \h 



c h 



B:=^F*{t)F*{t){hl\hl 

Y 



C:=Y,Fl(t)F Y (t)(h Y \h Y ). 



(20a) 



(20b) 



(20c) 



The log likelihood ratio of (18) can be analytically 
maximized over the amplitude parameters, resulting in 
the maximum likelihood ratio T -statistic [4], 



T : = lnA(x; s ML ) = -x^M^ x„, 



(21) 



where :— {M 1 }* J,V is the inverse of M^, 

i.e. M.^M.a^ — 6£, and takes the following form, 



B —C " ^ 

l_ I -C A 

D I B -C 

-C A 



(22) 



where D := AB — C 2 . It should be noted that the 
/"-statistic is the same as the square of coherent SNR 
(2F = Pcoh)i which has been previously used in litera- 
ture associated with coherent searches for inspiral signals 
with ground-based GW detectors [4, 5]. 
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IV. APPLICATION TO INSPIRAL SIGNALS 

So far our treatment of the ^-statistic could be equally 
applied to searching for GW signals from rotating neu- 
tron stars or inspiralling binaries of compact objects. Re- 
stricting ourselves to the case of non-spinning inspiral 
signals, the intrinsic parameters include {rj,M c }, where 
r] := ra\rail(m\ + rri2) 2 is the symmetric mass ratio 
and A4 C :— (mi + 7712)1^^ is the chirp mass. In ad- 
dition, there is an extrinsic parameter that can be effi- 
ciently maximized over but has not been in deriving the 
^-statistic, namely the coalescence time t c . This can 
be easily done with the use of the Fast Fourier Trans- 
form. The additional extrinsic parameters that must be 
searched over with separate filters are the sky locations 
{a,<5}, where a is the right ascension and 5 is the decli- 
nation. 

Although we restrict the derivation to the case of short- 
duration non-spinning inspiral signals, spins aligned with 
the angular momentum could easily be incorporated into 
the phase model and included as intrinsic parameters. 
This is because binaries in which the objects' spins are 
aligned with the angular momentum do not precess. 

In second generation GW detectors, the sensitive band 
of the detectors will start as low as 10Hz. A binary neu- 
tron star system's GW signal will be in the sensitive band 
of the detectors for < 17 minutes before coalescing, which 
amounts to a rotation of the Earth of < 0.07 radians. 
Thus, for a source's fixed sky location, the detectors can 
approximated as fixed djAt) « eK(t c ). With this approx- 
imation, the polarization weighted waveforms are given 
in the frequency domain as 



h\{f) 


= Fl 


hj(f) 


= Fl 


hj(f) 


= Fl 


hj(f) 





'h c (f), 



(23) 



The frequency domain intrinsic waveforms are given as 

h c (f) :=A(/)3?e^\ 
h s (f) :=A(/)3e^\ 



(24) 



where A(/) is the intrinsic amplitude of the waveform, 
^(f) is the phase of the waveform, and 5i and 3 denote 
operators that extract the real and imaginary parts, re- 
spectively. Each of these components has the following 
dependencies 



FX = F Y (t c ;a,S), 
F% = FZ(t c ;a,6), 
A(/) = A(/;M c ,r?), 



(25) 



NB: A(/) is typically defined to include ho and &(f) is 
typically defined (e.g., ['■)]) to include (f>o, however in this 



treatment, ho and <f>o are instead included as part of the 
amplitude parameters. The explicit expressions of A(/) 
and ty(f) according to the Stationary Phase Approxima- 
tion are expanded to Newtonian order in the amplitude 
and 3.5 PN order in the phase in Appendix A. 

The template waveforms for inspiral signals occupy a 
large bandwidth within the detectors, entering the sensi- 
tive band at the lower frequency cutoff f\ ow and extend- 
ing up to the frequency associated with the inner-most 
stable circular orbit /isco For these signals, the inner 
product between two waveforms is defined as 



/high 



Or 



:=43? 



flc 



(f)y Y *(f) 
s Y (f) 



df. 



(26) 



where /high is the upper cutoff frequency given by the 
smaller of the Nyquist frequency of that data or /iscOi 
x(f) denotes the Fourier transform of x(t), (.)* denotes 
the complex conjugate operator, and S Y (/) is the one- 
sided power spectral density (PSD) of detector Y. 



V. ^-STATISTIC METRIC DERIVATION 

The metric on the full set of parameters {A} includ- 
ing intrinsic and extrinsic parameters can be derived by 
expanding the log likelihood ratio (17) for a mismatched 
signal to second order in the parameter differences, AA, 

21n A(s(A); s(A + AA)) = 

(s(A)KA)) - (d a s(X)\d b s(X))AX a AX b 

+ G(AA 3 ), (27) 

where d a := d/d\ a is the partial derivative w.r.t. param- 
eter A a . In this notation, we will restrict the use of greek 
indices to the amplitude parameters (i.e., A p = AP) and 
for the metric subspace associated with the amplitude 
parameters. Using (27), we are led to the definition of 
the full metric g ab , which measures the fractional loss of 
IT , as 



9ab 



(d a s\d b s) 
(s\s) 



(28) 



Recalling that {d^h^dvA^hp) 
and using the signal model from (7), this metric can be 
decomposed into blocks 



1 



A^n v ^ A a h a0lJ A? 



(29) 



ah 



As stated before, in the above equation, the indices /i 
and v are associated with the amplitude parameter sub- 
space and the indices i and j are associated with the 
non-amplitude parameter subspace. The quantities TZ^i 



and h^ij are defined as 
V 



(dihn\djh v ) 



(30) 
(31) 
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The M-fw block is associated with derivatives of only 
the amplitude parameter subspace, the A a h a ^A^ block 
only with derivatives of the non-amplitude parameter 
subspace, and the A a lZ lia i block with derivatives of both 
subspaces. 

To obtain the metric for the ^-statistic, we can project 
out the dimensions associated with the amplitude sub- 
space [13, 21] 



{H,Ji,Ki,Gij}, and the notation x ■ y := ^yX f '\p' ', 
which denotes a sum over detectors,. The detectors' po- 
larization response vectors are defined as 



fpqY ._ pY p 

J ■ ± p ± c 



rpqY 



d,F Y F Y 



(37) 



(38) 



9ij = 9a - 9i a g afs gi3 r 



(32) 



Using the form of the full metric from (29), the F- 
statistic metric can be written as 



9ij A"M„,A'- ' 



(33) 



pqY 



d-F Y d F Y 



(39) 



where the derivatives of the detector polarization re- 
sponses diF Y x are given in Appendix B. Next, the de- 
tectors' waveform vectors {H,Ji,Ki,Gij} are defined 

as 



where the projected Fisher matrix O^vij is given by 

QpiHj = hiiuij — TZ afJ _. l M al3 lZf3v : j. (34) 

The two pieces of the projected Fisher matrix we re- 
fer to as the non-amplitude parameter subspace matrix, 
hfu/ij, and the amplitude subspace maximization correc- 
tion, 1Z ati iM a ^TZp V j. Similar to the derivation in Ap- 
pendix B of Ref. [13], after symmetrizing on (fJ.,1 1 ) and 



h/ivij takes the form 



(Ph 





4 

ij 



P?. 
Pi 

-h 







-pi 

\ lJ 
P 1 

p.3. 
'./ 



Pf. 

'j 



p?. 



(35) 



Although this looks identical to the derivation for ro- 
tating neutron star signals [13], one difference to keep in 
mind is that for inspiral signals there are additional terms 
hidden in these components. These additional terms are 
the result of the presence of an intrinsic parameter in the 
amplitude of the signal, what we refer to as the intrinsic 
amplitude. This can be seen in (25). The components of 



are given as 



-ft+-H, (36a) 



G Y := (h Y ddnA\h Y d ] \nA) + (h Y d l ^ Y (f)\h Y d j ^ Y (f)), 

(40) 

(41) 



H Y := (h Y \h Y ) 



Jf := (h Y \h Y dAnA) 



Kl := (h Y \h Y d^(f)), 



(42) 



(43) 



where the terms (h Y \h Y ), (h Y \h Y d l In A), 
{h Y \h Y d^ Y (f)), (/^lnAl^djlnA), and 

\h Y d^ Y {f)\h Y d^ Y (f)) are given in Appendix C. 
As referred to above, the additional terms for inspiral 
signals associated with derivatives of the intrinsic 
amplitude are contained in the G,j and J \ terms. 

Looking at the amplitude subspace maximization cor- 
rection, K ali iM. a ^Kp u j , IZ^i has the block form 



k, n t 

-Ki Hi 



where the blocks 1Z; and 7vL are defined as 



Ki := 



R} 1 

Rl 1 



Rf 



and Ki := 



R 13 

% 

R} 4 



Rf 

Rr 



(44) 



(45) 



P?i = f 



■ Gi 



+ 2 (ft> 



•J,+/ ? xx •J i + /?- x H, (36b) 



n + ) ■ Ji + i (/r + f* + ) ■ J i 



+ - (f +> 



fl + )-H, (36c) 



Hi = \ (ft X ft + )-Ki + \ (fj X ff + )'Ki. (36d) 



These components are defined in Appendix D. As noted 
in Appendix B of Ref. [13], Ki contains only terms with 
derivatives of the phase. However for the case of inspiral 
signals, IZi also contains terms with derivatives of both 
the antenna factors and the intrinsic amplitude. After 
using the symmetries of (22) and (44), and symmetrizing 
on (i,j), the final form for K all iM a ^TZp V j is 



TZ^M^K^j 



Qij Qij 



Qf 3 \ 



\Qt 



o -Q% Qh (>'■ 



o 



Q 



(46) 



ji V 



Above, we have introduced the detectors' polarization 
response vectors fujy the detectors' waveform vectors 



Explicit expressions for the Q components are given in 
Appendix D. 
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Combining the terms from (35) and (46), we find the 
projected Fisher matrix for inspiral signals has the same 
form as that of the low-frequency limit of rotating neu- 
tron star signals (i.e., Appendix B of Ref. [ ]), 



V 



v 



mij 
\ 





4 




\ 


4 

- m ij 






m v 


mj 







) 



(47) 



where mf- = P-j — Q\y Combining (33) and (47) gives 
the main result of this paper, namely, the coherent F- 
statistic metric for short-duration non-precessing inspiral 
signals. 

It should be noted that, as in the rotating neutron star 
case, although the F-statistic metric (33) has projected 
out the amplitude parameter subspace, it is still depen- 
dent on the amplitude parameters. This means that what 
has been derived is actually a family of metrics that de- 
pend on the the extrinsic parameters that enter the am- 
plitude parameters [13]. In order to produce a metric that 
is useful for choosing template points to cover the param- 
eter space, we must choose an averaging procedure. As 
an example, Prix takes the average of the eigenvalues of 
(A4 Q/3 t/ Q /3ijAA I AA J ) to produce an average metric. This 
is motivated by the fact that this matrix determines the 
extremal mismatches that can be obtained for any com- 
bination of amplitude parameters [13, 22]. 



VI. VERIFICATION 

With the F-statistic metric for short-duration non- 
precessing inspiral signals in hand, we can verify its 
performance by comparing the fractional loss of the F- 
statistic for mismatched signals to that predicted by the 
metric. We do this using a network of detectors corre- 
sponding to the locations and orientations of the LIGO 
Hanford, LIGO Livingston, and Virgo detectors. The 
PSDs we use for the LIGO detectors is the zero-detuning 
high-power advanced detector configuration [ ]. For 
Virgo, we use the advanced detector PSD [24]. The 
waveform model used for this work is the non-spinning 
restricted TaylorF2 PN approximation [25, 26], which 
is given in Appendix A. For computational reasons, we 
start the waveforms at a low-frequency cutoff of 40Hz, al- 
though our results should also be valid for other choices 
of the low-frequency cutoff. 

We perform our tests using the following intrinsic pa- 
rameters for the injected signal: mi = rri2 = 1AM@. The 
extrinsic parameters are: (a, 5) — (0,0), t c = 0, 0o = 0, 
if) = 0, cos l = 1, and V = 200 Mpc. The expected square 
coherent SNR for this signal is 2F = 9.8 2 . We check the 
metric by computing the match, both with and without 
maximization over time, while varying a single parame- 
ter. We do this for the two intrinsic parameters {M c , f?}, 
for the two extrinsic sky-location parameters {a,<5}, and 
also for the time parameter {t c }. Figure la shows how 



the match varies when the template's right ascension de- 
viates from the signal's value, shown as the vertical line. 
The metric reliably predicts the observed loss in F above 
~0.95. 

We are interested to see the effect that including 
derivatives of the detector responses has on the metric 
calculation. To do this, first we check the mismatches 
mfj from (47) associated with the F-statistic metric as 
a function of sky location, which can be seen in Fig. 2a. 
Figure 2b shows the portion of these mismatches that 
originates from the derivatives of the detector responses. 
We see that for the first three mismatches, this portion 
is typically an order of magnitude smaller than the full 
mismatch. As the fourth mismatch is already an order 
of magnitude smaller than the first three, including these 
terms is generally only a small correction to the metric. 
However, as we shall see, there are points in parameter 
space where this is not true. 

Finally, we check the effect of including the derivatives 
of the detector responses in the metric in an extreme ex- 
ample. We use the following intrinsic parameters for the 
injected signal: mi = ttt-2 = 1AM@. The extrinsic pa- 
rameters are: (a, S) = (0.785,-0.785), t c — 0, 0o = 0, 
ij) = 0, cos l = 0, V = 9.8 Mpc. The distance is an order 
of magnitude smaller than the previous comparisons in 
order to obtain an equal expected square coherent SNR 
for this signal, 2F = 9.8 2 . Figure 3 compares the pre- 
dictions from the metric derived with and without the 
derivatives of detector responses to the observed time- 
maximized fractional loss of F. We see that the predic- 
tions from the metric that includes the derivatives of the 
detector responses gives a substantially better match to 
the observed time-maximized fractional loss of F. How- 
ever, it should be noted that the detector network is much 
less sensitive to this point, which was chosen especially 
to show a large discrepancy between including versus not 
including those derivatives. For the majority of parame- 
ter space, the discrepancy is much smaller. 



VII. CONCLUSION 

In this work, we derive the coherent F-statistic met- 
ric associated with short-duration non-precessing inspiral 
signals. This metric, understandably, has very close ties 
to the coherent F-statistic metric associated with rotat- 
ing neutron star signals. However, in detail, there are 
several important differences. For one, inspiral signals 
have a larger bandwidth, hence the important single de- 
tector quantities are not the detectors' PSD values at 
a single frequency, but the integrated noise moments of 
the detectors' PSDs. Secondly, the signal model includes 
intrinsic parameters in the amplitude, which need to be 
properly accounted for in the metric derivation. 

Even though this derivation closely follows that for the 
rotating neutron star case, it includes previously ignored 
effects of the variation of the detector responses. If de- 
sired, this could easily be incorporated into the rotating 
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t c (ms) 
(e) 

FIG. 1: We show the fractional loss of J 7 as a function of parameter mismatches for the coherent analysis, using the 

Advanced LIGO and Advanced Virgo detectors, of an inspiral signal with component masses mi = mi = 1AM®, 
sky location (a, 8) = (0,0), and coalescence time t c = 0. The vertical line shows the true parameters of the injected 

signal. The solid lines show the observed fractional J 7 with and without maximization over time (where 
appropriate). The dashed lines shown the predicted fractional T from the coherent metric, without (gnu) & n d with 
(<7™^ xjr ) projection of the time dimension of the metric. The metric accurately predicts the fractional loss of T 
above a match of ^0.95. Panel (a), (b), (c), (d), and (e) show the mismatch associated variations of the right 
ascension, declination, chirp mass, symmetric mass ratio, and coalescence time, respectively. 
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FIG. 2: We show the separate mismatch components of the metric, mjy of (47), as functions of sky position in (a). 
The parts of associated with the derivatives of the detector response are shown in (b) . For the first three 
mismatches, these terms are typically more than an order of magnitude smaller than the full mismatches. 
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FIG. 3: We compare the time-maximized ^-statistic 
metric component with and without the derivatives of 
the detector responses, denoted as 3^ axjr and ggg ax:F ', 

respectively, for a specific set of parameter space 
coordinates: component masses mi = to 2 = IAMq, sky 
location (a, 5) — (2.7, 0.5), and coalescence time t c = 0, 

reference phase 4>q — 0, polarization angle ip = 0, 
inclination angle cost = 1, distance T> = 200 Mpc. The 
metric derived with the derivatives of detector responses 
better matches the observed time-maximized fractional 
loss of T . 



space (e.g., the sky space) that could be tiled separately. 
This would allow filters associated with different intrinsic 
parameters to be reused for the extrinsic parameters that 
still need to be searched in a tiled manner [11, 18]. 

Finally, because of the close ties between the metric 
and the projected Fisher matrix, it may be interesting 
to use the derivation here to determine the sky localiza- 
tion accuracy of a detector network, which could then be 
compared to the derivations of [14, 15, 27-29]. 
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neutron star coherent ^-statistic metric for a more com- 
plete picture of the sky-tiling problem. 

Important aspects that should be explored in the fu- 
ture include determining other ways that the amplitude- 
dependent metric, derived here, can be averaged [13] and 
applying the averaged metric to the template covering 
problem associated with coherent searches short-duration 
non-precessing inspiral signals. In order to efficiently per- 
form this search, it will need to be investigated how well 
the metric can be separated into an intrinsic parameter 
space (e.g., the mass space) and an extrinsic parameter 



Appendix A: TaylorF2 PN Waveform 

In this section we give the explicit formulae for the 
restricted SPA TaylorF2 inspiral waveform. As noted 
in Sec. IV, the inspiral waveform can be split into three 
pieces, a frequency-independent extrinsic amplitude (i.e., 
a function of only extrinsic parameters), a frequency- 
dependent intrinsic amplitude (i.e., a function of intrinsic 
parameters and frequency), and a phase piece that de- 
pends on intrinsic parameters, extrinsic parameters, and 
frequency. The extrinsic amplitude for a signal at dis- 
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tance T> is given by, 



h = 



5 1 
24 7r 2 / 3 r>' 



(Al) 



and the intrinsic amplitude for a signal with chirp mass 
M c is 



A(/) = M- 5 /3/-7/6. 



(A2) 



For convenience, we define A without the frequency de- 
pendence as 



A := AC 5/3 - 



(A3) 



The phase of the inspiral waveform, expanded to 3.5 PN 
order, can be written as 



V Y (f) = 2nft Y --+cb a 

+ I>/ ( " 5+J ' )/3 

6 

+ J2^Mf)f ( - 5+j)/3 , (A4) 

3=5 

where i y := t c — ry ■ fi/c is a time parameter that in- 
cludes the time of arrival of the end of the waveform at 
the geocenter t c and the sky-location-dependent correc- 
tion associated with a detector Y's location, and ipj and 
ijjj are the phase coefficients associated with the j/2 PN 
order. These phase coefficients are given by 



3M C 



5/3 



1287T 5 / 3 



V'2 



bM^rr 2 ^ /743 



V 

384tt 



84 



V"3 



-3^/ 3 X c - 2/ V 3/5 



(A5a) 
(A5b) 

(A5c) 



V>4 



3072tt 1 /3 

/ 3058673 5429 



V 7056 



7 



r? + 617/7' 



V, 5 = ^ f 7729 - 13, 
v 384 I 84 ' 



l + log( 6 3/2 7rX c 77- 3 / 5 , (A5e) 



Vj 384 V 84 



i>6 



TT^Ml^r)- 6 / 5 ( 11583231236531 
128 V 1564738560 ' 

,iS4S S + log [a^M]'^- 1 ' 
-3147553127 



640tH 



7 



451tt^ 



254016 



76055 o 127825 , 



576 



432 



06 



42 



(A5g) 
(A5h) 



^7 = 



5t: 5 / 3 M 2 c /3 v- 7 / 5 



V 

32256 

/ 15419335 75703 



V 336 



ry - 148097T 



(A5i) 



Any PN coefficients of 3.5 PN order or lower not defined 
above are identically zero. 



Appendix B: Explicit expressions for detector 
polarization response derivatives 

Based on the expressions for the detector polarization 
responses in Sec. II and [19], the derivatives of the de- 
tector polarization responses can be obtained in terms of 
derivatives of the polarization-independent basis tensors, 

'U- v ;>■:<■:. <>),■ (bi) 

These in turn can be written in terms of derivatives of 
the radiation frame basis vector, 

0«6? = {(0al5®f+f® (9a0 

-(daV)(E>rf-rf(g,(d a rf)y j , (B2a) 

d a 4 = {(5 |) Bjf+I® (d a rf) 

+(d a rf) <8 | + r/® (d a £)} lJ ■ (B2b) 

The explicit formulae for the derivatives of the radiation 
frame basis vectors with respect to the right ascension 
and declination are given as 



(A5d) 


d a 


h = (— cos 8 sin a, cos 8 cos a, 0), 


(B3a) 




dsn -- 


= (— sin 8 cos a, — sin 8 sin a, cos 8), 


(B3b) 






da£ = (cos a, sin a, 0), 


(B3c) 


(A5e) 




dsi= (0,0,0), 


(B3d) 






9 Q 77 = (sin 8 sin a, — sin 8 cos a, 0), 


(B3e) 


(A5f) 


dsfj = 


(— cos 8 cos a, — cos 8 sin a, — sin 8). 


(B3f) 
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Appendix C: Explicit expressions for inner product 
derivatives 



In this section we define the inner products that are 
needed for the coherent .F-statistic for short-duration 
non-precessing inspiral signals in terms of derivatives of 
intrinsic and extrinsic parameters and combinations of 
detector noise moment integrals, given by (El) of Ap- 
pendix E. 

The simplest inner product required, which contains 
no derivatives, is used by M.^ and H, 



(h Y \h Y ) = A 2 I(7,0y 



(CI) 



The inner product that contains a single derivative of the 
intrinsic amplitude and is used by J \ is 



(h Y \h Y dihiA) 



5 9 In M c A 2 T/n n ,y 



3 d\ i 



-A 1(7, 0) 



J 



(C2) 



The inner product that contains a single derivative of the 
phase and is used by Ki is 



(h Y \h Y d^ Y (f)) 



a 2 [%™m Y 



-E 



dX 



-1(12 -k,oy 



1(12 -fc,l) 



=:K 



(C3) 



Finally, the inner product that contains both two sin- 
gle derivatives of the intrinsic amplitude and two single 
derivatives of the phase is used by G !3 , 

(h Y dAnA\h Y d J \nA) + (h Y d^ Y (f)\h Y d^ Y (f)) =: G ly 
These individual inner products are given by 



Appendix D: Explicit terms associated with Q^vij 

The explicit formulae the amplitude subspace maxi- 
mization correction Q of the projected Fisher matrix, 
written in terms of the ic 4 components, are given as 



1 pl4^ 



DQl 3 = A(RfRf + R} 4 R 

C( | p21 nl 

— O JXj + H i tlj 



J3\fti it j T xtj tij 



pl3 pl4 
it, it, 



i^Ef) 



(Dla) 



DD 2 — A i P-^2 £?22 i p24 p24\ _i_ r>/ pl2 pl2 i pl4 pl4 
^ij — ^Hitj i% "r •"-» -"-j J "I" ^{ft-i ftj + itj 



- C(R\ 2 Rf 



RfRf 



R^Rf + RfRf), (Dlb) 



nn 3 — A ( E?21 p22 , p22 p21 , pl4 p24 , p24 pl4\ 
^ij — r -*S i% + itj itj -f- itj JXj ) 



11 pl2 

j 



pl2 pll , pl3 pl4 , pl4 pl3^ 



B(R] 1 R 

oil r>22 I r>22 tdII . r>12 r>21 , r>21 r>12 
- O [iti Itj ~r it,i JXj -+- -ft,; Jlj -+- ftj 

+ R?Rf + R 2i R? + 2R] i Rf), (Die) 



nr> 4 _ ^/p21 n24 i p24n21 pl4 p22 p22 pl4\ 
^ij ~ ftj + ftj ~ itj ftj ~ fti ftj ) 

I td( pll pl4 i pl4 pll Rl3pl2 pl2pl3~i 

+ tjlyft-i ftj T rtj ftj — fti ftj — itj -tlj j 

f~t( pll p24 , p24 pll , pl4 p21 , p21 pl4 
— O (JXi JXj + itj JXj + JX^ JXj + itj JXj 

pl3p22 p22 pl3 pl2pl4 pl4pl2\ mi,l\ 

— .ftj ttj — itj n,j — JX i JXj — rtj JXj ), [ULa) 



where we recall that A, B, C, D come from (22) for A4 
These Ri components are given by 



R] 1 = /++ H + /++ • Jj 



(D2a) 



(^^lnAl^^lnA) 



Ri 2 = fr-H+f^-j t 



(C5) 



and 



(/j y a I * y (/)|/ l y ^* y (/)) = A 2 



4ir 2 I(l,0y 



(dt Y d^ k d^ k dt Y \ y 

k y ' 
k v 7 

<9f/>fc ^ , , y 



y;^Mj(i7-fc-i,i) y 

^ 



k,l 



E 

fej 



1(17 -k -1,2) 



(C6) 



Rt = fr-H+r*Ju 



Rf = fr-H+r*.j u 



Rf = f ++ -K u 



R? = r*-Ki 



(D2b) 
(D2c) 
(D2d) 
(D2e) 
(D2f) 
(D2g) 



It is interesting to note that all of the Q components 
contain terms associated with derivatives of the detector 
polarization responses as well as terms associated with 
derivatives of the intrinsic amplitude. 
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Appendix E: Detector PSD Moment Integrals 

In this section define the noise moment integrals 
I(k, l) Y of detector Y's PSD, 



i(k,iy 



/high 



ln(/)/ 



-fc/3 



-df. (El) 



S Y (f) 



This is the same definition as in [9]. Based on 
the definition in (El), it is easy to see that pow- 
ers of the frequency in the inner-product can be 
manipulated as (h Y \n[(f)f- k ^ 3 \h Y \n l 2 (f)f- k2 / 3 ) = 
(h Y \h Y m *i+'2(y)/-(fci+fc 2 )/3) _ A 2 /(fc 1 + fc 2 + 7,Z 1 -H 2 ) y 
The detector PSD moments required for the metric cal- 
culation associated with restricted SPA TaylorF2 inspiral 
waveforms expanded to 3.5 PN order are 

(k,l)e {(1,0), (2,0), (3,0), (4,0), (5,0), (6,0), 

(7.0) , (8,0), (9,0), (10,0), (11,0), 
(12,0), (13,0), (14,0), (15,0), (17,0), 

(3.1) , (4,1), (5,1), (6,1), (7,1), 

(8.1) , (9,1), (10,1), (11,1), (12,1), 

(5. 2) , (6, 2), (7, 2)}. 
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Appendix F: Derivatives of PN Coefficients ^7^4 

Here we give explicit expressions for the derivatives of 
the PN coefficients associated with the phase in terms of 
the symmetric mass ratio r\ and the chirp mass M, as 
computed in [9]. First, the derivatives with respect to 
M, d vA 



dMtpo = 
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d n tp 7 = 



3225677 12 / 5 

/ 15419335 
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V 48 



7570377 + 44427?7' 



(F16) 



Any derivatives of the coefficients of 3.5 PN order or lower 
not defined above are identically zero. 

The phase also contains a term associated with the 
relative time shift between the arrival of a signal at a 
detector's location compared to the arrival of the signal 



at the geocenter. These derivatives are as follows 
Y r Y ■ d a h 



(F17a) 



ry ■ dsn 



(F17b) 



where d a n and dsn are given by (B3a) and (B3b), re- 
spectively. 
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